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Abstract 

We analyze the rate of detection of minimal supersymmetric neutralino dark matter 
in germanium, sapphire and sodium iodide detectors, imposing cosmological and recent 
accelerator bounds including those from b — S7 decay. We find, in contrast with 
several other recent analyses, that although the b — S7 constraint reduces the number 
of viable models, models still remain where the counting rate in solid state detectors 
exceeds 10 kg~^ day~^. 

The recent observation by the CLEO collaboration of the b — ?■ S7 decay [1] has stirred 
interest in the possible bounds obtainable for supersymmetric models that contribute 
to this process [2]. Some authors [3, 4] have gone one step further and analyzed the 
consequences of this and other accelerator bounds for the predicted rates in experiments 
aimed at detecting neutralino dark matter. One problem with even the so-called "minimal" 
supersymmetric standard model (MSSM) is that it contains a large number of parameters 
which in principle are free, although several of them are constrained by various naturalness 
arguments and experimental bounds on flavor changing interactions etc. In order to reduce 



the number of parameters, some assumptions related to supersymmetric grand unification 
of coupling constants and masses are usually made. Even more restricted models have 
been suggested (e.g., [5]), where the simplest possible structure is chosen at the GUT 
scale and the electroweak symmetry breaking is assumed to be achieved only through 
radiative corrections to the Higgs potential when the GUT parameters are run to lower 
energies using the renormalization group equations. 

As expected, making the model more restrictive gives more restrictive bounds on the 
experimental signals of the lightest supersymmetric particles, both in accelerator experi- 
ments and in dark matter detection experiments (if supersymmetry is related at all to the 
dark matter problem of the galactic halo). In particular, it was shown in [3, 4, 5] that 
rates for both direct and indirect detection of neutralino dark matter generally become 
very small in this type of models. At this point it should be realized, however, that ev- 
idence is mounting that whatever the effective low-energy supersymmetric model is, it is 
unlikely to be of such a simplified form. For instance, in models derived from superstrings 
it is not clear if there is an intermediate GUT scale at all. Even if there is, threshold 
corrections at the GUT scale may be important (for a recent discussion of some of these 
problems, see [6]). In fact, recent analyses of the MSSM phenomenology have appeared 
in which the universality condition of the scalar masses at the GUT scale is relaxed, es- 
pecially for the Higgs sector [7]. It has also recently been suggested [8] that to reproduce 
the measured value of a, at low energies, it may even be necessary to modify one of most 
frequently used relations, that of the gaugino mass parameters (we follow the notational 
convention of [9]): 

Ml = I tan^ e^M2 ^ O.SMj, 

M2 = ■ ?r M3 ^ O.3M3. ^ > 

The upshot of this is that although the simplistic models have the virtue of being 
predictive, and due to small rates largely irrelevant for present-day dark matter detectors, 
the price paid is the lack of generality. 

If one adopts a phenomenological approach and allows for a more general variation 
of parameters in the MSSM, still of course consistent with experimental bounds and giv- 
ing correct low-energy symmetry breaking, one may be tempted to define a probability 
measure in parameter space. That is, one may want to translate the number of allowed 
models giving, e.g., a range of detection rates to a "probability". This method is often 
used implicitly as one typically scans a large number of models and presents results as 
density of points in some parameter plane. It may then be tempting to assign values where 
there is a large density of model points a higher probability than those with small density. 
We will argue, illustrated by some simple examples, that this is not meaningful to do, 
since different parametrizations of the models can give completely different results. 

The conservative approach we propose is to regard the whole range of outcomes of a 
calculation as a priori equally probable, irrespective of the parametrization. This means 
that really only upper and lower limits can be given. Another consequence of this philos- 
ophy is that it becomes important to scan the model parameter space accurately enough 
to discover also the very extreme outcomes of the calculation, since they will be used for 
setting the bounds. 
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We started this study with the intention of relaxing the various GUT or radiative 
breaking conditions to explore the full range of possibilities in a generic minimal super- 
symmetric model. We found, to our surprise, that even without relaxing those assumptions 
there still were models in which the dark matter detection rates in solid state detectors 
were quite high. We think that this discrepancy with published results may be due to 
insufficient scanning of parameters space in the earlier works. 



The minimal supersymmetric standard models 

We work in the framework of the minimal N = 1 supersymmetric extension of the standard 
model defined by, besides the particle content and gauge couplings required by supersym- 
metry, the superpotential 

W = eij (-e^Y^ii^/ - d^Yzpqi^/ + u^jYc,qi^| - /x^l^|) (2) 

and the soft supersymmetry-breaking potential 

Vsoit = eij (-e^A^Y^li^/ - d^AzpY^pqi^/ + u^Ac,Yc,qi^| - S/x^l^| + h.c.) 
+Hi*mlHi + HiTmlHi, 

+^MiBB + (W^W^ + 2W+W-) + ^Msgg. (3) 

Here i and j are SU(2) indices (ei2 = +1), Y's, A's and M's are 3x3 matrices in 
generation space, and the other boldface letter are vectors in generation space. 

The one-loop effective potential for the Higgs fields in the dimensional reduction renor- 
malization scheme then follows (see [10]): 

V = hIh{*H{ + nlHiTHi, - B/j, (eijHlHi + h.c.) 

+\{g' + g'') [Hi*Hi - W*Hl)' + \g' |^r^||' + AFi , (4) 

with = /x^ + , = /x^ + m\ and the one-loop corrections 



(5) 



The supertrace is defined as ^tvf{M'^) = Ei C'i(-l)^*'(2si + ^)f{mj) where d is the 
color degrees of freedom and Si is the spin of the i^^ particle. 

Electroweak symmetry breaking is caused by both Hi and acquiring vacuum ex- 
pectation values, 

{Hl) = vu {Hl) = V2, (6) 

with g'^{vl + = 2m^, with the further assumption that vacuum expectation values of 
all other scalar fields (in particular, squark and sleptons) vanish. This avoids color and/or 
charge breaking vacua. 
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At the tree level, the correct symmetry breaking pattern results if 



fil + fil>2\Bfi\, (7) 

to guarantee a potential bounded from below, and 

filfil < B^ii\ (8) 

to have a saddle point and not a minimum at vi = V2 = 0. At the one-loop level, the 
potential is automatically bounded from below thanks to the logarithms in the one-loop 
contributions, and the origin is not an extremum. 

The minimization conditions of the potential (4) can be written as 



= ^/2v 
= V2v 



'nl cos2 15 -nl sin2 /3 + + g'^y cos 2/3] + ATi , 
+ sin 2/3 - S/i] + AT2 , 

where v"^ = v\ + tan/3 = v^/vi and ATi, AT2 are one-loop tadpole contributions [10]. 
These minimization conditions allow one to trade two of the Higgs potential parameters 
/xf, /x| and Bfj, for the Z boson mass rri^ = 1(5^ + fl''^)(^i + ^2) ^"^^ ^^e ratio of vevs 
tan/3. The third parameter can further be reexpressed in terms of the mass of one of the 
physical Higgs bosons. 

When diagonalizing the mass matrix for the scalar Higgs fields, besides a charged and 
a neutral would-be Goldstone bosons which become the longitudinal polarizations of the 
and Z gauge bosons, one finds a neutral CP-odd Higgs boson A, two neutral CP-even 
Higgs bosons Hi^2 and a charged Higgs boson H^, which will play an important role in 
the calculation. Choosing as independent parameter the mass of the CP-odd Higgs 
boson, the masses of the other Higgs bosons are given by 



nasses of the other Higgs bosons are given by 

cos^ /3 + m| sin^ /3 + AA^^i - sin/3 cos/3(m^ + m|) + AA^^2 A 
- sin/3cos/3(m^ + m|) + AA^Ii sin^ /3 + m| cos^ /3 + AA^^j / 



mjj± = m\ + + A±. (11) 

The quantities A.Ai^j and A± are the one-loop radiative corrections coming from virtual 
(s)top and (s)bottom loops, calculated within the effective potential approach as in [11]. 
Diagonalization oi Ai^ gives the two CP-even Higgs boson masses, 2 > and their mixing 
angle a {—n/2 < a < 0). 

The neutralinos Xi are linear combination of the neutral gauge bosons B, W3 and of 
the neutral higgsinos H°, H^. In this basis, their mass matrix 



is diagonalized analytically to give four neutral Majorana states, 

X° = ZiiB + Zi2W^ + ZizH[ + Zi4^2°, (13) 

the lightest of which, to be called %, is then the candidate for the particle making up the 
dark matter in the universe. 

The charginos are linear combinations of the charged gauge bosons and of the 
charged higgsinos , H2 ■ Their mass mass terms are given by 



{W- H-) U2+J^^"'" ^^^^ 



Their mass matrix. 



M,.= [fl "?). (15) 



is diagonalized by the following linear combinations 

Xi = UiiW- + Ui2H-, (16) 
Xt = ViiW+ + Vi2Ht. (17) 

We choose det(C/) = 1 and U*M.y.±V^ = diag(m~±, m~±) with non-negative chargino 
masses m~± > 0. 

When discussing the squark mass matrix including mixing, it is convenient to choose 
a basis where the squarks are rotated in the same way as the corresponding quarks in the 
standard model. We follow the conventions of the particle data group [12] and put the 
mixing in the left-handed c?-quark fields, so that the definition of the Cabibbo-Kobayashi- 
Maskawa matrix is K = Vi Vj, where Vi (V2) rotates the interaction left-handed u-quark 
(c?-quark) fields to mass eigenstates. For sleptons we choose an analogous basis, but due 
to the masslessness of neutrinos no analog of the CKM matrix appears. 

We then obtain the general 6 X 6 u- and c?-squark mass matrices: 

^ / + mt m„ + Dlj^l -mt (A^, + /x* cot/3) \ 
" \ -{Au + Hcot(3)mu + m„mt + D^jjl ; ' ^ ^ 

^ ( KtM^K + lAimt + Di^l -mt (At, + ^^* tan/3) \ 
V -(Azp+/itan/3)md Mf, + + D^^jl / ' 

and the general sneutrino and charged slepton masses 

Ml = Ml + Dl^l (20) 

2^ fMi + memt+D£^l -mt(A^ + /x*tan/3) \ .^l) 
" V -(A£ + /itan/3)me M| + mtnie + 7 ' ^ ' 



Here 



D{i = m|cos2/3(T3/-e/sin2^^), (22) 
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Drr = ml COS 2l3efsm'^e^. (23) 

In the chosen basis, iriu = diag (rrzu, rric, mt), m<i = diag (mj, nis, rrzb) and me = diag(me,mjti, ttIt). 

The slepton and squark mass eigenstates fk (pk with k = 1,2,3 and ej., % and c?fe 
with k = 1, . . ., 6) diagonalize the previous mass matrices and are related to the current 
sfermion eigenstates and via (a = 1, 2, 3) 

ha = j2hn^£. (24) 

k=i 

fRa = j^hn^ (25) 

k=l 

The squark and charged slepton mixing matrices Tul,r, ^dl,r and Tel,r have dimension 
6x3, while the sneutrino mixing matrix Ti,l has dimension 3x3. 

For simplicity, and to get in touch with published papers [3], we then make a simple 
ansatz for the up-to-now arbitrary soft supersymmetry-breaking parameters: 

Au = diag(0, 0, At) 
Ad = diag(0, 0, Ab) 
Ae = 

Mq = Mu = Md = Me = Ml = mol. 

This allows the squark mass matrices to be diagonalized analytically. For example, for the 
top squark one has, in terms of top squark mixing angle 0^, 

rM = ^¥r = cos r|i = -r*}^ = sin (27) 

Notice that the ansatz (26) implies the absence of tree-level flavor changing neutral 
currents in all sectors of the model. It is not however the more general ansatz for the 
absence of tree-level FCNC's, which would only demand that the trilinear couplings A^, 
the soft mass matrices M? and also K^MgK be diagonal, and not necessarily equal to each 
other. Notice also that this is typically not what is obtained in low-energy supergravity 
models with a universal scalar mass at the grand-unification (or Planck) scale, in which 
the running of the scalar masses down to the electroweak scale generates off-diagonal terms 
and tree-level FCNC's in the squark sector. 



(26) 



The b ^ S7 branching ratio 

In the calculation of the branching ratio for b — ?■ S7 we follow [13]. In the standard 
model, the rate is given by 

r(b ^ ST) = ^j^, (28) 



with 



K^^K.^Xtwfiixtw), (29) 
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where e and g are evaluated at the b mass scale, and the function /i arises from the loop 
integration: 

8a;3-3a;2_i2a; + 7 + 6a;(2-3a;)log(a;) , , 

^^(") = 36(. - 1)4 • 

Here and in the following i/nfij. 

For a general supersymmetric model, contributions from charged Higgs bosons H~ , 
charginos gluinos g and neutralinos x° have to be added to A^. In our case, only 
the standard model contribution Aw and those for H~ and x~ enter, since the gluino 
and neutralino contributions vanish because of our ansatz for the squark mass matrices 
(absence of tree-level FCNC's). 

Ajj- and A^- can be written as (see, e.g. [13]) 

2 

^H- = ^/f 2 ^*s^tb XtH \cot' /3 /i (xtH) + /2 (xtH)] , (31) 



2 6 

EE 

c=l fe=l 



m.- 

Xc 



(32) 



where ~ ~± , and On~ ~± , are the left- and right-handed down-quark-chargino-up-squark 

■^^feXc ^'^kXc ^ 

vertex constants given by 



9lu,xU = E {VcI^urY^ - aV:,^^^) Kid 

i=l 



The additional loop integration functions entering are 

5a;2 - 8a; + 3 + 2(2 - 3a;) log(a;) 



/2W 



and 



/3(^) 



6(a;-l)3 
-7x^ + 12a; - 5 - 2a;(2 - 3a;) log(a;) 



(33) 
(34) 

(35) 
(36) 



6{x - 1)3 

It is important for our result that the loop integration functions fi{x) [i = 1,2,3) are 
decreasing positive functions of x [14]. 

The b — 7- S7 branching ratio is finally evaluated as 

r(b ^ S7) 



with 



BR(b^s7) = BR(b^ce.)^^^^^^^^ 
BR(b cev) = 0.107 



(37) 
(38) 
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and 



r(b ^ cev) = -^lifebl V(a;c6)A(m6, x,k), (39) 



where 

p{x) = 1 - 8a;2 + Sa;^ - a;* - 24a;^ log(a;) (40) 
is a phase space factor, and 

2 

A(m, a;) = 1 - — a,(m)/s(a;) (41) 

OTT 

is a QCD correction factor for the semileptonic decay {fs{xcb) — 2.4 [15]). 

A complication in the analysis is given by possibly large QCD corrections to the b — ?■ 
S7 decay, which at present are plagued by theoretical uncertainties (see e.g. ref. [16]). 
Therefore we preferentially opted for the tree-level value described in this section. In the 
discussion we will comment on the effect of including QCD corrections according to the 
procedure outlined in ref. [13]. 

The direct detection rate 

The rate for direct detection of galactic neutralinos, integrated over deposited energy 
assuming no energy threshold, is 

R = J2Nin^{<^i^v), (42) 

i 

where Ni is the number of nuclei of species i in the detector, is the local galactic 
neutralino number density, CTj^ is the neutralino-nucleus elastic cross section, and the 
angular brackets denote an average over v, the neutralino speed relative to the detector. 

We take the local galactic neutralino velocity distribution as a truncated gaussian, 
which in the detector frame moving at speed vq relative to the galactic halo reads 



A/cut uvQcr 



{u - voY 



2(72 



exp 



2^^ 



(43) 



for Uesc < < \/^esc + (^o + ^cut)^ and zero otherwise, with u = \/v^~-\-v^ and 



Numerically, we have taken the halo line-of-sight velocity dispersion a =120 km/s, the 
galactic escape speed v^ut = 600 km/s, the relative Earth-halo speed vq = 264 km/s 
(a yearly average) and the Earth escape speed Vesc = 11.9 km/s. We have adopted a 
local dark matter density m^n^ = 0.3 GeV/cm^ whenever the calculated neutralino relic 
density n^^/i^ > HgaLDM^i^) the density of dark matter in galactic halos averaged over the 
entire universe. This is meant to represent the minimum value for which neutralinos could 
make up the totality of the galactic dark matter. When n^^/i^ < HgaLDM^i^) we have scaled 
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proportionally to ^x/^galDM- The value to choose for ^gaLDuh"^ is quite uncertain, 
both because of uncertainties in the density and extension of galactic halos and because 
of the poorly known relation between the universally-averaged and the local dark matter 
densities. Just for comparison, we have chosen the same value ^gaLDuh"^ = 0.025 as in 
ref. [3], but we call to the attention of the reader that values smaller by one order of 
magnitude would still be acceptable. 

The neutralino-nucleus elastic cross section can be written as 

where rriiy. is the neutralino-nucleus reduced mass, q is the momentum transfer and Gi^i^q^) 
is the effective neutralino-nucleus vertex. In this calculation we write 

GUl') = AlFl{q')Gl + 4AM{<l')Gl, (46) 

(often A? appears as X^J[J + 1) ) and assume gaussian nuclear form factors [17] 

Fs{q') = FA{q') = exp{-q' /6fi'), (47) 

Ri = (0.3 + 0.89AV^)fm, (48) 

which should provide us with a good approximation of the integrated detection rate [18], 
in which we are only interested. Since the non-zero-spin nuclei we consider have an un- 
paired proton, we consider the scalar and axial neutralino-proton vertices Gs and Ga, 
a neutron admixture giving a contribution that can be neglected. Using heavy-squark 
effective lagrangians [19], we get 

Gs= (^9) I E ^^^^f - 1 ^ gLg.xim.x, j (49) 



g=u,d,s,c,b,t \h=Hi,H2 ^ k=l Qk 



and 



^" \ 8fi — ^ — ^ ^ 

g=u,d,s \ ^ k=l Ik / 

The gi's are elementary vertices involving the particles indicated by the indices, and they 
read 

^ f (^^^2 - ffy^xi) (--^x3Cosa + Z^4sina) , for 
^^^^ \ {gZx2- 9yZxi){Zx3sma + Z^icosa), for if 2, ^ ' 

_ / -YqCosa/^/2, ioT Hi, 
^'^^^ - 1 +F, sin a/A for ^2, ^ ' 



9Zxx 



' - Zl,) (63) 

9LqkXq = 9LL^QL+9RL^QR, (55) 

9Rqkxq = 9lrTql + 9RR^QR, (56) 



2 COS 6w 
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with 



9LL = -^(TsqgZ^2 + ^gyZ^i^ , (57) 

QRR = V2eqgyZ^i, (58) 

/ -YqZ^s, forg = u, c, t, , , 

9LR = 9RL = \ V ry t 11 (59) 

I -YqZ^^, for g = d,s,b, 



and 



= J ->2, forg = u,c,t, 
I niq/vi, lor g = d, s, b. 

Numerically, we have taken [20] 

mu(uu) = 0.023mp, md(dd) = 0.034mp, (61) 
ms{ss) = 0.14mp, mc(cc) = mb(bb) = mt(tt) = 0.595mp, (62) 

and [21] 

Au = 0.77, Ad =-0.49, As = -0.15. (63) 
Moreover, we have used 

Ail = 0.35, A^^ = 0.041 and A| = 0.007, (64) 

according to the odd-group model [22]. 

One should be aware that both our choice of nuclear form factors and of neutralino- 
nucleon vertices and the numerical values adopted for the nucleon matrix elements are 
at best approximate. A more sophisticated treatment (see discussion and references in 
[23]), would however, beyond still presenting theoretical and calculational weaknesses, 
change the rate values by much less than the spread due to the unknown supersymmetric 
parameters.^ 

Sampling of supersymmetric parameter space 

When solving the minimal supersymmetric standard models defined above we let the 
universal (weak-scale) sfermion mass parameter mo vary between 100 and 3000 GeV, 
the soft supersymmetry breaking parameters Af, and At between — 3mo and 3mo, tan/3 
between 1.2 and 50, the pseudoscalar mass ttia between its experimentally allowed lower 
bound and 1000 GeV, and we fix the top quark mass at rrit = 175 GeV. 

For the model scan, two alternative sets of parameters are used: (1) the usual /x and 
M2, both varied logarithmically in the interval [—5000, 5000] GeV, and (2) the mass of the 
lightest neutralino and its gaugino fraction Zg = Z^^ + Z^2i varied logarithmically in 



^The b ^> S7 formulas in ref. [23] and in the accompanying computer code neutdriver seem to contain 
some errors. After correction and fixing of some bugs, neutdriver gave results in qualitative agreement 
with ours. We thank G. Jungman for providing us with the source code. 
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[-5000, 5000] GeV and linearly in [0.00001, 0.99999]. To start with, we enforce the gaugino 
mass unification conditions (1). 

With this definition of the MSSM models (which we note is the same as in [3]), we 
generate model parameters randomly within the bounds mentioned above and proceed to 
calculate the mass spectrum, the couplings, the Z-boson invisible width, etc. We keep 
only models that satisfy the accelerator constraints given in the 1995 Review of Particle 
Properties [12]. In addition, we drop models that violate the 95% C.L. limits from the 
CLEO experiment [1] 1.0 • 10"^ < BR(b sj) < 3.4 • 10"^. We ask for 4,000 models 
satisfying the accelerator constraints including b — ?■ S7. 

Then for each model allowed by the accelerator constraints we calculate the relic den- 
sity of neutralinos Q^-h"^. We use the formalism in ref. [24] to carefully treat resonant 
annihilations and threshold effects, keeping finite widths of unstable particles, including 
all two-body annihilation channels of neutralinos. The annihilation cross sections used 
were derived using a novel helicity projection technique [25], and were checked against 
published results for several of the subprocesses. Only models that would not overdose 
the universe, i.e. in which n^^/i^ < 1, are considered cosmologically viable. 

In fig. 1 we show the /X-M2 location of the viable models that survive all experimental 
bounds, including the n^^/i^ and b — ?■ S7 constraints. Fig. 1(a) refers to the /X-M2 sampling, 
and fig. 1(b) to the m^-Zg sampling. 

At this point we compute integrated direct detection rates for pure germanium (^^Ge), 
sapphire (AI2O3) and sodium iodide (Nal) detectors, which are representative of devices 
presently under research or development. The rates resulting in experimentally and cos- 
mologically viable models are presented in fig. 2 for a Ge detector. Again, fig. 2(a) refers 
to the /X-M2 sampling, and fig. 2(b) to the m^-Zg sampling. 

The huge spread of possible rates, even at fixed neutralino mass, is evident. This 
hinders the predictability of the class of MSSM models we consider, and is one of the 
driving reasons for introducing more restrictive, and so more predictive, relations among 
the model parameters. 

A fact should however be noticed: in the two samplings, no points have interesting 
detection rates, say above 1 event/kg/day in Ge. Can we conclude that in the class of 
models we consider the b — ?■ S7 constraint is so strong to exclude interesting detection 
rates? 

Compare figs. 2(a) and 2(b): the aesthetic appearance of the clouds of points is quite 
different. This is only due to the different a priori probabilities used in the two samplings. 
It is obvious that Nature, if she has chosen supersymmetry, has realized just one of the 
models appearing as dots in our figures. By the same token, we have no reason to prefer 
an a priori probability measure over another. And so we must not be lead astray by 
attaching a probability to the points in the figures. We must not turn sentences like 'most 
of the models give low detection rates' into statements on the relative likelihood of high 
and low detection rates. The same applies of course to histograms derived from such 
samplings. And if comparison of figs. 2(a) and 2(b) is not convincing enough, it should 
become inescapable in the following. 

We perform two special scans, of 500 models each, demanding to consider only models 
in which the scalar piece of the neutralino-proton cross section is larger than a tenth of 
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the corresponding Dirac neutrino-proton cross section. This in order to pick out points 
with the highest detection rates. In the first scan we sample in /X-M2 as before but restrict 
the range of the Higgs pseudoscalar mass to ttia G [0, 60] GeV. In the second scan, we 
sample in m^-Zg with the restricted ttia range and further demand G [800, 1200] GeV 
and Zg G [0.01, 0.99]. The results of these special scans for Ge, AI2O3 and Nal are shown 
in fig. 3, combined with (the top parts of) those of the previous samplings. Remarkably, 
the high-rate zones, empty before, are now filled with points. Particularly striking is the 
concentration of points around — 1000 GeV, which obviously comes from the second 
special sampling. In fig. 4 we also show the dependence of the counting rates on the 
gaugino fraction Zg for the combined sample. The upper band, corresponding to the two 
special scans, is almost flat. The slightly higher rates at Zg > 0.5 are essentially due to a 
weaker or even absent H-suppression of the galactic neutralino density. 

We remind that all points shown in figs. 3 and 4 are compatible with the experimental 
and cosmological constraints mentioned above, which include the b — ?■ S7 bounds. And 
so we conclude that, as far as only these constraints are considered, there are viable 
models with an integrated counting rate as large as 10 events/kg/day in Ge, and even 
100 events/kg/day in Nal. These models might well be already probed (and excluded) by 
current dark matter searches. 



Discussion 

As we have shown, contrary to the results of [3], we find models with as high Ge 
counting rate as 10 events/kg/day which still do not violate the b — ?■ S7 bound. 

In fig. 5 we plot the b — ?■ S7 rate as a function of the mass. The Ge detection 
rate decreases from the upper left to the lower right panel. The vertical band of points 
at relatively low mjj± belongs to the two special samples, and shows best at high Ge 
rates. The curve at relatively large mjj±, most pronounced at very low Ge rates, comes 
from the naive /X-M2 and m^-Zg samplings combined, and describes the naively-expected 
relationship between detection rate and charged Higgs boson mass. The two groups of 
points are detached, but this is only an artifact of the sampling procedure. 

The reason for having an acceptable b — ?■ S7 from a relatively light (100 - 200 GeV) 
happens to be a cancellation between the charged Higgs contribution and the contribution 
from the lighter of the two charginos [14]. The latter may be sizable at large tan/3 when 
the top-squark mixing is substantial and the lightest chargino is dominantly a charged 
higgsino xt ~ relatively light (IM2I 3> \n\ > mw)- Under these conditions, the 

dominant chargino contribution is well approximated by 



-sign(/x%) 



eg 



167r2m^ 



mw 



^/2\m~- 1 cos /3 
X2 



V ~-f,JV~\f,+ sin^rcos^r /sfaJr — /sfaJr 



(65) 



This charged higgsino contribution can be negative and effectively cancel the W boson and 
charged Higgs boson contributions, which are always positive, when the top squark mixing 
angle 0^ and the parameter /x have the same sign. In terms of the soft-supersymmetry 
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breaking parameters this amounts to 

niAt + /xcot/3) < 0, (66) 

or in common instances in which At = O(mo) 3> /x/ tan/3 to roughly fiAt < 0. In our class 
of models it is possible to satisfy the previous condition (66) for both /x positive and /x 
negative, because we are free to choose the sign and magnitude of At. This freedom is lost 
in models imposing additional theoretical constraints, for example in no-scale models or 
in models with a flat Kahler manifold [A = OoTA = B — m respectively at the unification 
scale). 

We have also performed an analysis including QCD corrections to the b — ?■ S7 am- 
plitudes according to the prescription in ref. [13]. For a given model, the BR(b — ?■ S7) 
is generally larger than the tree level value, so some of the models that were previous 
viable have too large a QCD-corrected BR(b — ?■ S7). However, other models that had 
too low a tree-level b — ?■ S7 decay rate become viable when the QCD corrections increase 
BR(b —7- S7). So our conclusions on the existence of models with high counting rates and 
acceptable BR(b — ?■ S7) remain valid even after including QCD corrections. 

We already mentioned that we have chosen ^gaLDuh"^ = 0.025 just for the sake of 
comparison with ref. [3], and that the actual value to use might even one order of magnitude 
larger or smaller. Were it smaller, many models in which the detection rate is suppressed 
just because of a too small relic density would add to the number of models with important 
detection rates. Fig. 6 shows the predicted rates in Ge versus the calculated neutralino 
relic density Q^-h"^, for the combined sample of models. The two naive samples fill the 
triangular shape at the bottom, and the two special samples are the band and cloud in the 
upper parts. Models to the right of the solid vertical line overdose the universe (and have 
been plotted to illustrate the trend of R versus Q^h"^). Rates to the left of the vertical 
dashed line [D,^ = rigai.DM) are H-suppressed by our simple-minded prescription for the 
local galactic neutralino density at low D,^. Lowering Hgai.DM shifts the tip of the 'volcano' 
upwards and towards the left so that right-hand side continues straight to the left, and 
raises the cloud of specially-sampled models to still higher rates. This would also lower the 
typical neutralino mass (here 1000 GeV) for which calculated detection rates are highest. 
Notice that in some of the interesting models, the neutralino relic density is larger than 
what is needed for them to fill up galactic halos, i.e. D,^ > rigai.DM- 

Some of the models shown in the figures may indeed be already excluded by current 
negative searches of halo neutralinos. We have not tagged these models as excluded 
because a serious analysis would require a more detailed calculation of the predicted rates, 
involving different differential rates for each experiment, sophisticated nuclear form factors, 
quenching factors, etc. Such an analysis is out of the scope of the present paper. 

We set out to go beyond the restrictive supergravity models and examine neutralino 
detection rates in a general minimal supersymmetric model. Relaxing our ansatz (26) on 
the sfermion masses would demand the consideration of several phenomenological bounds 
from flavor changing neutral currents, as e.g. the K°K° and B°B° mass differences, the 
electric dipole moment of the neutron, etc. We have started by relaxing the request of a 
universal gaugino mass, eqs. (1), inspired by a recent suggestion to reproduce the measured 
value of as at low energies [8]. We have replaced eqs. (1) by Mi — O.3M2 — M3 and run 
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through our calculation again. The plots we obtain do not differ qualitatively from those 
we presented, and the quantitative differences are slight. For these reasons, we do not show 
them here. Our conclusions on the existence of supersymmetric models with acceptable 
b — 7- S7 branching ratios and important detection rates remain the same. 

Conclusions 

In this paper we have shown that even with the same kind of supersymmetric models 
as in [3], and enforcing the constraints from b — ?■ S7, we obtain a range of predictions for 
direct dark matter search experiments that are much less restrictive. In particular we have 
found models in which the integrated counting rates in Ge detectors are calculated to be 
higher than 10 events/kg/day, and even higher in Nal detectors. The basic reason behind 
the compatibility of these models with the b — ?■ S7 constraints is the non-existence of a 
lower bound on the charged Higgs boson from b — ?■ S7 measurements, since in the presence 
of a large top squark mixing the W and contributions to the b — ?■ S7 amplitude may 
be effectively canceled at large tan/3 by the contribution from a light charged higgsino. 
This cancellation may occur at both positive and negative values of /x. 

We have argued that with the present ignorance of the origin of an effective low-energy 
supersymmetric theory one should allow the phenomenological parameters to vary over 
as large a range as possible without imposing unnecessary relations between parameters 
(but still restricting them according to various experimental bounds coming from particle 
physics and cosmology). Although we have kept some simplifying assumptions, the large 
range of the predicted rates is remarkable. As dark matter detectors improve, we expect 
this range to be successively narrowed. 

We have stressed that the density of points in the plots depends on the assumed a 
priori distribution of the parameters. One should not be misled in thinking that rate 
values in zones where there are more points are more probable than those in which there 
are few. To illustrate this, we have shown the rates obtained by trading /x and M2 with 
the neutralino mass and gaugino fraction Zg. The aesthetic appearance of the plots 
is indeed different. To further stress our point, we have also presented the rates for 
an a priori distribution which privileges light charginos and light charged Higgs bosons. 
With this choice, high detection rates with acceptable b — ?■ S7 branching ratios look 
"generic." It is therefore apparent that no probability should be attached to the plotted 
point distributions (or to histograms derived from them), and that the figures can only 
illustrate possible neutralino detection rates. With this caveat, the detection rates are 
highest for a neutralino mass around 1000 GeV, preferentially more gaugino than higgsino, 
but these neutralino characteristics are sensitive to the prescription for the H-suppression 
of the galactic neutralino density. 

Naturally, one should consider properties that are invariant under an arbitrary change 
of the a priori parameter probabilities, like the maximum (or minimum) values of the 
quantities of interest (the detection rates in our case). Unfortunately, a thorough and 
fine scanning of parameter space is computationally very expensive and an alternative 
analytical extremization seems prohibitively complicated. We therefore have to leave the 
following question open: are there in fact additional, allowed points in the empty regions 
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of our plots? 
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Figure captions 

Figure 1: Location of experimentally and cosmologically viable models in the /X-M2 
plane for (a) the /X-M2 sample and (b) the m^-Zg sample. 

Figure 2: Integrated direct detection rate R off ^^Ge versus neutralino mass for (a) 
the fi-M2 sample and (b) the m^-Zg sample. 

Figure 3: Integrated direct detection rate R in ^^Ge, AI2O3 and Nal detectors versus 
neutralino mass for the combined sample. 

Figure 4: Integrated direct detection rate R in ^^Ge, AI2O3 and Nal detectors versus 
neutralino gaugino fraction Zg for the combined sample. 

Figure 5: Branching ratio for b — ?■ S7 decays versus charged Higgs boson mass mjj± 
for the combined sample. 

Figure 6: Integrated direct detection rate R off ^^Ge versus neutralino relic density 
n^^/i^ for the combined sample. 
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